Robust genetic analysis of the X-linked anophthalmic (Ie) mouse

Anophthalmia (missing eye) describes a failure of early embryonic ocular development. Mutations in a relatively small set of genes account for 75% of bilateral anophthalmia cases, yet 25% of families currently are left without a molecular diagnosis. Here, we report our experimental work that aimed to uncover the developmental and genetic basis of the anophthalmia characterising the X-linked Ie (eye-ear reduction) X-ray-induced allele in mouse that was first identified in 1947. Histological analysis of the embryonic phenotype showed failure of normal eye development after the optic vesicle stage with particularly severe malformation of the ventral retina. Linkage analysis mapped this mutation to a ~6 Mb region on the X chromosome. Short- and long-read whole-genome sequencing (WGS) of affected and unaffected male littermates confirmed the Ie linkage but identified no plausible causative variants or structural rearrangements. These analyses did reduce the critical candidate interval and revealed evidence of multiple variants within the ancestral DNA, although none were found that altered coding sequences or that were unique to Ie. To investigate early embryonic events at a genetic level, we then generated mouse ES cells derived from male Ie embryos and wild type littermates. RNA-seq and accessible chromatin sequencing (ATAC-seq) data generated from cultured optic vesicle organoids did not reveal any large differences in gene expression or accessibility of putative cis-regulatory elements between Ie and wild type. However, an unbiased TF-footprinting analysis of accessible chromatin regions did provide evidence of a genome-wide reduction in binding of transcription factors associated with ventral eye development in Ie, and evidence of an increase in binding of the Zic-family of transcription factors, including Zic3, which is located within the Ie-refined critical interval. We conclude that the refined Ie critical region at chrX: 56,145,000–58,385,000 contains multiple genetic variants that may be linked to altered cis regulation but does not contain a convincing causative mutation. Changes in the binding of key transcription factors to chromatin causing altered gene expression during development, possibly through a subtle mis-regulation of Zic3, presents a plausible cause for the anophthalmia phenotype observed in Ie, but further work is required to determine the precise causative allele and its genetic mechanism

Monoallelic variants resulting in substitutions of MAB21L1 Arg51 Cause Aniridia and microphthalmia

Classical aniridia is a congenital and progressive panocular disorder almost exclusively caused by heterozygous loss-of-function variants at the PAX6 locus. We report nine individuals from five families with severe aniridia and/or microphthalmia (with no detectable PAX6 mutation) with ultrarare monoallelic missense variants altering the Arg51 codon of MAB21L1. These mutations occurred de novo in 3/5 families, with the remaining families being compatible with autosomal dominant inheritance. Mice engineered to carry the p.Arg51Leu change showed a highly-penetrant optic disc anomaly in heterozygous animals with severe microphthalmia in homozygotes. Substitutions of the same codon (Arg51) in MAB21L2, a close homolog of MAB21L1, cause severe ocular and skeletal malformations in humans and mice. The predicted nucleotidyltransferase function of MAB21L1 could not be demonstrated using purified protein with a variety of nucleotide substrates and oligonucleotide activators. Induced expression of GFP-tagged wildtype and mutant MAB21L1 in human cells caused only modest transcriptional changes. Mass spectrometry of immunoprecipitated protein revealed that both mutant and wildtype MAB21L1 associate with transcription factors that are known regulators of PAX6 (MEIS1, MEIS2 and PBX1) and with poly(A) RNA binding proteins. Arg51 substitutions reduce the association of wild-type MAB21L1 with TBL1XR1, a component of the NCoR complex. We found limited evidence for mutation-specific interactions with MSI2/Musashi-2, an RNA-binding proteins with effects on many different developmental pathways. Given that biallelic loss-of-function variants in MAB21L1 result in a milder eye phenotype we suggest that Arg51-altering monoallelic variants most plausibly perturb eye development via a gain-of-function mechanism

Identification and functional modelling of plausibly causative cis-regulatory variants in a highly-selected cohort with X-linked intellectual disability.

Identifying causative variants in cis-regulatory elements (CRE) in neurodevelopmental disorders has proven challenging. We have used in vivo functional analyses to categorize rigorously filtered CRE variants in a clinical cohort that is plausibly enriched for causative CRE mutations: 48 unrelated males with a family history consistent with X-linked intellectual disability (XLID) in whom no detectable cause could be identified in the coding regions of the X chromosome (chrX). Targeted sequencing of all chrX CRE identified six rare variants in five affected individuals that altered conserved bases in CRE targeting known XLID genes and segregated appropriately in families. Two of these variants, FMR1CRE and TENM1CRE, showed consistent site- and stage-specific differences of enhancer function in the developing zebrafish brain using dual-color fluorescent reporter assay. Mouse models were created for both variants. In male mice Fmr1CRE induced alterations in neurodevelopmental Fmr1 expression, olfactory behavior and neurophysiological indicators of FMRP function. The absence of another likely causative variant on whole genome sequencing further supported FMR1CRE as the likely basis of the XLID in this family. Tenm1CRE mice showed no phenotypic anomalies. Following the release of gnomAD 2.1, reanalysis showed that TENM1CRE exceeded the maximum plausible population frequency of a XLID causative allele. Assigning causative status to any ultra-rare CRE variant remains problematic and requires disease-relevant in vivo functional data from multiple sources. The sequential and bespoke nature of such analyses renders them time-consuming and challenging to scale for routine clinical use

FRA2A is a CGG repeat expansion associated with silencing of AFF3

Folate-sensitive fragile sites (FSFS) are a rare cytogenetically visible subset of dynamic mutations. Of the eight molecularly characterized FSFS, four are associated with intellectual disability (ID). Cytogenetic expression results from CGG tri-nucleotide-repeat expansion mutation associated with local CpG hypermethylation and transcriptional silencing. The best studied is the FRAXA site in the FMR1 gene, where large expansions cause fragile X syndrome, the most common inherited ID syndrome. Here we studied three families with FRA2A expression at 2q11 associated with a wide spectrum of neurodevelopmental phenotypes. We identified a polymorphic CGG repeat in a conserved, brain-active alternative promoter of the AFF3 gene, an autosomal homolog of the X-linked AFF2/FMR2 gene: Expansion of the AFF2 CGG repeat causes FRAXE ID. We found that FRA2A-expressing individuals have mosaic expansions of the AFF3 CGG repeat in the range of several hundred repeat units. Moreover, bisulfite sequencing and pyrosequencing both suggest AFF3 promoter hypermethylation. cSNP-analysis demonstrates monoallelic expression of the AFF3 gene in FRA2A carriers thus predicting that FRA2A expression results in functional haploinsufficiency for AFF3 at least in a subset of tissues. By whole-mount in situ hybridization the mouse AFF3 ortholog shows strong regional expression in the developing brain, somites and limb buds in 9.5-12.5dpc mouse embryos. Our data suggest that there may be an association between FRA2A and a delay in the acquisition of motor and language skills in the families studied here. However, additional cases are required to firmly establish a causal relationship

Psip1/p52 regulates posterior Hoxa genes through activation of lncRNA Hottip

Long noncoding RNAs (lncRNAs) have been implicated in various biological functions including the regulation of gene expression, however, the functionality of lncRNAs is not clearly understood and conflicting conclusions have often been reached when comparing different methods to investigate them. Moreover, little is known about the upstream regulation of lncRNAs. Here we show that the short isoform (p52) of a transcriptional co-activator—PC4 and SF2 interacting protein (Psip1), which is known to be involved in linking transcription to RNA processing, specifically regulates the expression of the lncRNA Hottip–located at the 5’ end of the Hoxa locus. Using both knockdown and knockout approaches we show that Hottip expression is required for activation of the 5’ Hoxa genes (Hoxa13 and Hoxa10/11) and for retaining Mll1 at the 5’ end of Hoxa. Moreover, we demonstrate that artificially inducing Hottip expression is sufficient to activate the 5’ Hoxa genes and that Hottip RNA binds to the 5’ end of Hoxa. By engineering premature transcription termination, we show that it is the Hottip lncRNA molecule itself, not just Hottip transcription that is required to maintains active expression of posterior Hox genes. Our data show a direct role for a lncRNA molecule in regulating the expression of developmentally-regulated mRNA genes in cis

Clinical and molecular consequences of disease-associated de novo mutations in SATB2

Purpose: To characterize features associated with de novo mutations affecting SATB2 function in individuals ascertained on the basis of intellectual disability. Methods: Twenty previously unreported individuals with 19 different SATB2 mutations (11 loss-of-function and 8 missense variants) were studied. Fibroblasts were used to measure mutant protein production. Subcellular localization and mobility of wild-type and mutant SATB2 were assessed using fluorescently tagged protein. Results: Recurrent clinical features included neurodevelopmental impairment (19/19), absent/near absent speech (16/19), normal somatic growth (17/19), cleft palate (9/19), drooling (12/19), and dental anomalies (8/19). Six of eight missense variants clustered in the first CUT domain. Sibling recurrence due to gonadal mosaicism was seen in one family. A nonsense mutation in the last exon resulted in production of a truncated protein retaining all three DNA-binding domains. SATB2 nuclear mobility was mutation-dependent; p.Arg389Cys in CUT1 increased mobility and both p.Gly515Ser in CUT2 and p.Gln566Lys between CUT2 and HOX reduced mobility. The clinical features in individuals with missense variants were indistinguishable from those with loss of function. Conclusion: SATB2 haploinsufficiency is a common cause of syndromic intellectual disability. When mutant SATB2 protein is produced, the protein appears functionally inactive with a disrupted pattern of chromatin or matrix association

Identification and Validation of a Putative Polycomb Responsive Element in the Human Genome

<div><p>Epigenetic cellular memory mechanisms that involve polycomb and trithorax group of proteins are well conserved across metazoans. The <i>cis</i>-acting elements interacting with these proteins, however, are poorly understood in mammals. In a directed search we identified a potential polycomb responsive element with 25 repeats of YY1 binding motifthatwe designate PRE-PIK3C2B as it occurs in the first intron of human <i>PIK3C2B</i> gene. It down regulates reporter gene expression in HEK cells and the repression is dependent on polycomb group of proteins (PcG). We demonstrate that PRE-PIK3C2B interacts directly with YY1 <i>in vitro</i> and recruits PRC2 complex <i>in vivo</i>. The localization of PcG proteins including YY1 to PRE-PIK3C2B in HEK cells is decreased on knock-down of either <i>YY1</i> or <i>SUZ12</i>. Endogenous PRE-PIK3C2B shows bivalent marking having H3K27me3 and H3K4me3 for repressed and active state respectively. In transgenic <i>Drosophila</i>, PRE-PIK3C2B down regulates mini-white expression, exhibits variegation and pairing sensitive silencing (PSS), which has not been previously demonstrated for mammalian PRE. Taken together, our results strongly suggest that PRE-PIK3C2B functions as a site of interaction for polycomb proteins.</p></div

Interaction of PcG proteins with PRE-PIK3C2B in HEK cells.

<p>The experiments were carried out in HEK293 ells. A-Diagrammatic representation of the endogenous PRE-PIK3C2B region with primers (arrows) used for PCR and the 25 mer repeating unit (filled box). B and C- ChIP with antibodies indicated. Input is 20% of sonicated chromatin. D-Quantitative PCR for PRE-PIK3C2B following ChIP in control cells and cells transfected with siRNAYY1 and siRNA SUZ12. E-Interaction of EZH2 with PRE-PIK3C2B estimated by qPCR in absence and presence of YY1siRNA: F-Localization of EED and the H3K4me3 at PRE-PIK3C2B in absence and presence of siRNAYY1. Error bars: S.E.M of assay in triplicate, * p<0.05, **p-value<0.005,n = 3. The qPCR profile for positive and negative controls are shown in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0067217#pone.0067217.s003" target="_blank">Figs. S3</a> and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0067217#pone.0067217.s004" target="_blank">S4</a>.</p

Effect of Pho mutation on PRE-PIK3C2B mediated repression.

<p>A- Expression of miniwhite mRNA by qPCR. The genetic background is indicated below the histograms.* p<0.05, **p<0.005. B- comparison of eye pigmentation between transgenic flies with and without PHO mutation. ΔPI indicates the lines where PRE-PIK3C2B is flipped out.</p